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Abstract 

We present an improved determination of the up- and down-quark distributions in the proton us¬ 
ing reeent data on eharged lepton asymmetries from W- gauge-boson produetion at the LHC and 
Tevatron. The analysis is performed in the framework of a global fit of parton distribution fune- 
tions. The fit results are eonsistent with a non-zero iso-spin asymmetry of the sea, x(d-u), at 
small values of Bjorken x ~ 10“"^ indieating a delayed onset of the Regge asymptoties of a vanish¬ 
ing (J - M)-asymmetry at small-x. We compare with up- and down-quark distributions available in 
the literature and provide accurate predictions for the production of single top-quarks at the LHC, 
a process which can serve as a standard candle for the light quark flavor content of the proton. 


I. INTRODUCTION 


The hadro-production of the electroweak gauge bosons W- and Z is an important reaetion, 
whieh ean be measured with very high preeision in hadron eollider experiments, sueh as Tevatron 
or the Large Hadron Collider (LHC). The available theoretieal predictions match the precision of 
the experimental data thanks to the known radiative corrections to next-to-next-to leading order 
(NNLO) in QCD and to next-to leading order (NLO) in the electroweak sector of the Standard 
Model. The process of W-- and Z-boson production has, therefore, become known as a so-called 
standard candle process, because it provides valuable constraints on the parton content of the 
colliding protons. In a global fit of parton distribution functions (PDFs) precision data on W-- and 
Z-boson production allow to separate the quark flavors, in particular up- and down quarks, and to 
disentangle the sea and valence quark PDFs. 

Recently, new data sets have become available from the LHC [ii and the Tevatron experiment 
D0i,i based on measurements of the rapidity distributions of charged leptons produced in the 
decays of the -bosons or charged-lepton pairs from the Z-boson decays. The kinematic reach 
to large rapidities makes the data analysis very interesting for small-x physics, particularly in 
the context of the LHC. In the past, the proton structure at small x has been studied in detail 
at the HERA collider. However, the potential for quark flavor separation has always been poor 
due to the underlying deep-inelastic scattering (DIS) process probing dominantly flavor singlet 
states. Successful improvements in the determination of the strange quark PDF have recently 
been undertaken dH] with the help of new neutrino-nucleon DIS data, but the non-strange sea 
at small x could not be separated with existing data thus far. The latter, therefore, is routinely 
described by a parameterization with Regge-like asymptotics in global PDF fits. The new DO data 
on electron and muon asymmetries [S 0 as well as the LHCb data [jH] on the measurement of 
the forward VF“-boson cross-section in pp collisions can help to resolve this ambiguity due to 
the particular kinematics in the forward region. In addition, the CMS measurement of the muon 
charge asymmetry [HI] assists in the improved determination of light parton PDFs due to very high 
accuracy of the data. 

In this letter we study the impact of the new data on charged lepton asymmetries from W-- 
boson production and their rapidity distributions on the global fit of PDFs by ABM dfij-Isl]. Taking 
the ABM 12 fit [[sl] based on the world DIS data, measurements of Drell-Yan (DY) dimuon pro¬ 
duction from fixed targets and early LHC data on W-- and Z-boson production as a reference the 
impact of individual new data sets is quantified and the resulting shifts in the up- and down quark 
PDFs are documented. The fit results turn out to be consistent with a non-zero iso-spin asymmetry 
of the sea, I(x) = x{d-u), at small values of Bjorken x ~ 10“"^. This indicates a delayed onset 
of the Regge asymptotics of a vanishing asymmetry I{x) at small x. Particular care is taken to 
control potential correlations of the J-quark distribution at large x and the J-quark distribution at 
small X. The results are compared with earlier determinations of the asymmetric sea x{d - u) at 
small X obtained in the global fits CTIO [@], CT14 JT^, MMHT iflUl or the so-called unbiased fit 
NN3.0 [121 as well as with the recent HERAFitter analysis [ 131. 

The fit results of the present analysis for the u- and J-quark distributions are used to provide 
accurate predictions for the production of single top-quarks at the LHC. This process furnishes a 
very sensitive test for the light quark flavor content of the proton. The respective cross sections 
are almost complete to NNLO in QCD for t-channel single top-quark production iHl and mea¬ 
surements at the LHC have been conducted at = 7 TeV by ATLAS ifl^ and at = 8 TeV 
by CMS ifT^ . With the new PDFs of the present fit we obtain updated predictions for the single 
(anti)top-quark cross sections CTt and ai as well as their ratio Rt = ctiIctj. The latter quantity is par- 
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ticularly sensitive to PDFs as uncertainties largely cancel in the ratio. As a result, we demonstrate 
in this letter that single top-quark production can serve as a standard candle process at the LHC, 
cf. [El. 


II. SET-UP OF THE ANALYSIS 


The present analysis is based on the update of the ABM 12 global fit [[80 including recent data 
on the neutrino-induced dimuon production H and with the addition of recent data on W-- and 
Z-boson production collected by the LHC and Tevatron experiments. This new input allows to 
improve the determination of the quark distributions, and, in particular, the separation of the u- 
and J-quarks. At first instance, these are the LHCb H and DO [0] data on forward VT“-boson 
production providing a constraint on the sea quark iso-spin asymmetry x{d - fi) in a range of x, 
which has not been probed in earlier experiments. The LHCb statistics of Ref. [H is collected in 
the muon channel, with the muon pseudo-rapidity in the range of 2 4.5, which is sensitive 
to parton momentum fractions x down to x ~ 10“^. The same applies to the new DO data on the 
electron-positron asymmetry [0] covering the electron/positron pseudo-rapidity ? 7 ^in the range of 
-3.2 3.2. The upper margin of this range is smaller than one of the LHCb data yO, however due 
to the difference in the collision energy of Tevatron and the LHC, the DO data probe the quark 
distributions down to x ~ 10“"^ as well as the LHCb ones. The forward kinematics of the LHCb 
and DO experiments is also sensitive to the quark distributions at large x and, in particular, can be 
used to pin down the J-quark distribution, which is poorly known in this region. The final DO data 
on the muon charge asymmetry [0] span the region of |? 7 ^| <2. Therefore they put a constraint on 
the quark iso-spin asymmetry mostly at x > 0.001. The updated CMS data on the muon charge 
asymmetry [El] with < 2.4 are even more limited in the parton momentum coverage. However, 
they provide important experimental input due to unprecedented accuracy. We also add to the 
analysis recent data on the Z-boson production in the forward region ilsL llOll obtained by the 
LHCb experiment at the collision energies of 7 and 8 TeV. These data provide a complementary 
constraint on the PDF combination and allow further improvements in disentangling the quark 
species at small and large x. On the other hand, earlier data from CMS on the electron charge 
asymmetry and from LHCb on forward -boson production llElll are removed from the fit in 
view of their low statistical significance as compared to the more recent data of these experiments. 
The deuteron data by the NMC, BCDMS, and SLAC experiments employed in the ABM12 fit are 
also not included into the present analysis in order not to be exposed to uncertainties related to 
modeling of the nuclear effects. 

The collider data on W-- and Z-boson production used in the present analysis are summarized 
in Table [I| This set is essentially different from the respective data used in the ABM 12 fit. Indeed, 
many high-statistical samples are added now, while only the ATLAS data on W-- and Z-boson 
production [l22n and the LHCb data on Z-boson production in the electron channel ll23ll are used 
in both cases. 

In order to take advantage of this input to constrain the PDFs fully exclusive QCD calculations 
of the VF“-boson production process are required, because the data sets of Tableware obtained in 
a restricted phase space due to limited detector acceptance and the /Z-boson event selection 
criteria with a cut on the lepton transverse momentum imposed. Fully differential QCD pre¬ 
dictions for lepton rapidity distributions are implemented up to NNLO in the publicly available 
code FEWZ (version 3.1) [|2412511 . FEWZ (version 3.1) provides a convenient capability to estimate 
uncertainties in the cross sections due to the PDFs. The data sets newly included into our analysis 
are compared with the NNLO FEWZ predictions based on the ABM12, CTIO, CT14, MMHT, and 
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Experiment 

ATLAS 

CMS 

DO 

LHCb 

(TeV) 

7 

7 

1.96 

7 

8 

Einal states 

1T+ ^ Z+v 

W- tv 

z tr 

1T+ ^ 77+y 

IT- ^ 77-y 

1T+ ^ 77+y 

IT- ^ 77-y 

tT+ ^ e+y 

IT- ^ e-y 

1T+ ^ 7/+y 

IT- ^ //-y 

Z ^ e+g- 

Reference 

\m 

ill 

[3J 

[4] 

im 

09] 

Cut on the lepton Pj 

P^j > 20 GeV 

> 25 GeV 

/^ > 25 GeV 

> 25 GeV 

P^ > 25 GeV 

P^ > 20 GeV 

Luminosity (1/fb) 

0.035 

4.7 

7.3 

9.7 

1. 

2. 

NDP 

30 

11 

10 

13 

31 

17 


TABLE I: The different data samples of VT*- and Z-boson production from LHC and Tevatron used in the 
present analysis. 


NN3.0 PDFs in Figs. [B H [3] and H 

The predictions agree in general with the data ^. However, in places random fluctuations be¬ 
yond the uncertainties appear. In particular, those fluctuations are observed for the CMS and LHCb 
data on the lepton asymmetry, cf. Fig.[3l This prevents to obtain an ideal value of for these two 
samples. In the LHCb sample one specific data point at 7]^ = 3.375 deviates from the general trend 
most significantly. It is worth noting that the shape of final state radiation (FSR) corrections to 
account for QED corrections in W-- and Z-boson decays applied to the data of Ref. S demon¬ 
strates an anomalous irregularity precisely at this value of 77 ^, cf. Fig. [51 In contrast, in the earlier 
LHCb analysis of the muon asymmetry Kill] the FSR corrections employed exhibit a smooth shape 
and for such a shape the fluctuation in the pseudo-rapidity distribution at rj^ = 3.375 vanishes. The 
treatment of FSR in LHCb differs between the respective analysis. In Ref. [0] the LSR corrections 
are computed with the PHOTOS Monte Carlo while Ref. [18] has determined LSR corrections 
from the mean of simulations done with Herwig-i-i- [l27n and PythiaS KSil . Having no possibility 
to validate this issue, we simply drop the problematic points at 77 ^ = 3.375 in the LHCb VF“-boson 
data sample from the fit. This reduces the value of by some 10 units and brings it down to a 
reasonable level, although the LSR corrections still may need clarification for the overall LHCb 
kinematics of Ref. ifl^ and other publications of this experiment. In particular, this concerns the 
uncertainty in the LSR correction, which has been recently estimated as +0.5%, which is much 
bigger than the value of ±0.2% quoted in the earlier LHCb analysis of Ref. 

In line with this revision we also expand the LSR correction uncertainty in the LHCb Z-boson 
production data of Ref. ll23n to the same value of +0.5%. However, the early LHCb data of 
Ref. [|23n are anyway in a certain tension with more recent LHCb Z-boson production samples 
collected in different decay channels and at different collision energies. Indeed, the electron data 
of Ref. ll23n collected at the center of mass (c.m.s.) energy ^fs = 1 TeV overshoot (undershoot) 
the muon ones of Ref. ifl^ at the same ^[s when the pseudo-rapidity is bigger (smaller) than 
77 ~ 3, while the trend of LHCb electron data iflsll collected at = 8 TeV is similar to the one of 


' The LHCb data of Ref. |Q] on VL'^-boson production undershoot the predictions at 77 ^ ~ 2. However they also come 
in agreement once the trigger efficiency is fixed illSll . cf. Fig. [2] 


3 























Ref. UJ], cf. Fig. [I 


Experiment 

EHCb 

v^(TeV) 

7 

8 

Einal state 

Z ^ e+e- 

Z ^ 

Z ^ e^e- 

Reference 

[23J 

\m 

[19] 

NDP 

9 

17 

17 

x\ 

9.4 

- 

- 

xl 

- 

11.3 

- 

x] 

- 

- 

19.1 

xl 

16.4 

14.8 

- 

4 

16.2 

- 

21.4 

xl 

- 

10.5 

18.3 


TABLE II: The values obtained in the present analysis for the LHCb data on Z-boson production lIlSL 


ISLllJ], where index i = 1 ... 6 denotes variants of the fit based on different combinations of those data. 


This tension is quantified by the value of obtained in variants of our analysis, whieh in- 
elude different eombinations of the LHCb Z-boson data without adding the other data sets listed 
in Table U in order to provide a elean eomparison. A good value of;^f^ is only aehieved when one 
single LHCb Z-boson data set is ineluded into the fit (variants #1,2,3 of Table UH). A slightly in- 
ereased value appears only for the eleetron sample at = 8 TeV demonstrating fluetuations 
beyond the uneorrelated errors (of. Fig. [I]). However, upon oombination of the eleetron sample 
at = 7 TeV with the muon one dX ^fs = 1 TeV and the eleetron one at = 8 TeV (variants 
#4,5) the value of;^f^ inoreases substantially. For the oombination of muon data dX ^fs = 1 TeV and 
the eleetron ones at = 8 TeV samples (variant #6) on the other hand the value of remains 
essentially the same as for the fits with single samples. As a oonsequenoe of this observation we 
exclude the LHCb electron sample it^ at Vi = 7 TeV from our analysis. Moreover, the errors in 
those data are also substantially bigger than in the ones of Refs. HEi. 

The theoretical predictions displayed in Figs. [HO[3] and |4] demonstrate some random fluctua¬ 
tions due to the numerical integration of the NNLO QCD matrix elements in FEWZ. These fluctua¬ 
tions are inessential for the purpose of qualitative comparison. However, they cause an additional 
contribution to the value of x^ when quantifying the agreement of those theory predictions with 
the data. To suppress this contribution to a marginal level one has to provide an integration accu¬ 
racy better than the data uncertainties, i.e., at the level of several %o for the high-statistics W-- 
and Z-boson data samples used in our analysis. To achieve such an accuracy one has to generate 
2-10^ FEWZ histories for each bin in pseudo-rapidity distributions. 

The results obtained in this way are illustrated in Fig. 0 on the example of the DO data on 
electron-positron asymmetry [0]. The integration accuracy displayed in Fig. 0 is smaller than the 
experimental uncertainties, particularly in the forward region, where the data are less accurate. At 
this rate the value of x^ is only marginally sensitive to the integration uncertainties in the QCD 
matrix elements and the generation of additional FEWZ histories is not necessary. However, the 
typical CPU time necessary to accomplish these computations is quite substantial, of the order 
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(9(1 O^h) for each data set of Table HI Therefore they cannot be repeated iteratively, with variations 
of the PDFs in the fit. Instead, we compute the prediction only once, but for each eigenvector 
encoding the uncertainties in an initial PDF set. Afterwards, upon varying the PDFs in a fit, 
we obtain the W-- and Z-boson cross sections corresponding to the current PDF values as an 
interpolation between the entries of this eigenvector prediction grid (cf. Ref. Jsl] for more details). 
Conceptually, this approach is based on the algorithm employed in the code fastNLO for the 
fast computation of differential distributions of jets in hadro-production with account of the NLO 
QCD corrections [|2^ . The main difference with fastNLO is the selection of the PDF basis used 
to generate the prediction grid. In our case it is much smaller than the one of fastNLO. However it 
is sufficient once the PDF variations due to new data included are within their uncertainties. As an 
advantage, the computing time necessary for the grid generation is greatly reduced as compared 
to the case of a general PDF basis employed in fastNLO. 

To meet the increased statistical significance of the new data sets used, cf. Table H we allow 
more flexibility in the fitted PDFs by releasing a previously imposed constraint on the sea quark 
iso-spin asymmetry 

I(x) = x(ds -Us) ~ , ( 1 ) 

where Sud ~ 0.5 due to Regge model arguments. In the present analysis the boundary conditions 
for the u- and d sea quark distributions at the scale of /i^ are chosen 

xus{x,nl) = Us{x,nl) = Aus(l - , ( 2 ) 

Xds{x,nl) = ds{x,ldl) = Ads(l - ^ ( 3 ) 

respectively, where 

Pusix) = il+y-ius^nx)il+yi,usx), ( 4 ) 

Pdsix) = (1 +y-i 4 s\nx){\ +y\ 4 sx). (5) 

In this representation the parameters aus,ads and y-i,us,7-i4s control the low-x PDF behavior. 
Since they are fitted independently, as well as the normalization parameters Aus,Ads, the small-x 
shape of the sea quark iso-spin asymmetry in the parameterization Eqs. Q, ([3]) is unconstrained. 
The remaining PDF shapes for the valence, strange, and gluon distributions used in the present 
analysis are the same as the ones of the ABM 12 fit. 

Having checked the statistical significance of the data sets in Table U we find that the fitted 
value of the parameter y-i4s is comparable to 0 within uncertainties. Furthermore, the value of;^f^ 
remains essentially the same when it is fixed at 0. Therefore, we impose the constraint y-i 4 s = 0 
for all variants of the fit considered. The sea quark iso-spin asymmetry at the initial scale of 
yuo = 3 GeV obtained in this way is displayed in Fig. |7]m comparison with the earlier ABM 12 
determination. Due to the relaxed functional form in Eqs. @, ([3]) its central value goes lower 
than the ABM 12 one at x < 0.01. This range of kinematics is mostly controlled by the forward 
VF*-boson data of DO [0] and EHCb [|2], which clearly prefer a negative value of the asymmetry at 
X ~ 10“^. Meanwhile the central W-- and Z-boson data obtained at the EHC [Hi 112] demonstrate 
the same trend, although with less statistical significance (cf. Eig. |7]). In this context it is worth 
noting that the description of the ATEAS data in the present analysis is improved as compared to 
the ABM 12 one (cf. Table ITVl in Section |ni|), evidently due to the relaxed shape of the sea quark 
distributions at small x. At x < 10“^ the shape of /(x) demonstrates turnover. Therefore, the Regge 
asymptotics still may recover at even smaller values of x, say x ~ 10“^. However, this observation 
is based on an unsophisticated extrapolation of Eqs. (|2|), (|3]) to those kinematics, which is not 
controlled by existing data. 


5 


The forward VT“-boson data also control quark distributions at large x and, in particular, allow 
to constrain the d/u ratio in the absence of deuteron-target DIS data employed for this purpose in 
the ABM 12 fit. The shape of the d/u ratio determined in the present analysis is in agreement with 
the ABM 12 result (cf. Fig.|7]). This confirms the validity of the Kulagin-Petti approach OOll used 
for modeling nuclear effects in the ABM 12 fit. A general form of the boundary conditions for the 
valence u- and J-quarks defining the d/u ratio at large x is the same as for the sea quarks and reads 


xuyix,iul) = Auv(l - x )^"^, ( 6 ) 

xdv{x,iJ.l) = Advil - , (7) 

respectively. However, in this case the polynomials Puvix) and Pdvix) read 

Puvi^) — ^uv'^yi,uv^'^y2,uv^ +T3,t<v-^ 5 (8) 


and 


Pdvi^) — ^dvyi,dv^y2,dv^ > 


(9) 


respectively, allowing more flexibility for the large-x PDF shape. The coefficient y 34 v could not 
be constrained by the data used in the ABM 12 fit, therefore it was set to 0. The fitted value of 
y 3 ,dv obtained using the data of Table His also comparable with 0 within uncertainties. So y 3 ,dv 
is set to 0 in the present analysis as well. The potential of the collider data in the determination 
of the d/u ratio at large x is comparable to the one of the deuteron data used in the ABM 12 fit. 
The DO electron data demonstrate the best sensitivity to dju preferring a somewhat lower value as 
compared with the LHCb ones, cf. Fig. |71 It is worth noting that the combination of the proton- 
target DIS and Z-boson production data also provide some constraint on the d/u ratio at large x 
(cf. the shaded area in the right panel of Fig. |7]). Such a determination is consistent with the results 
of our combined analysis, although being much less accurate. 


Experiment 

DO 

LHCb 

(TeV) 

1.96 

7 

Einal state 

1T+ ^ e+v 

1T+ ^ e+y 


1 

t 

1 

> 

1 

T 

1 

Reference 

[4] 

[18] 

NDP 

13 

14 

x\ 

9.5 

- 


- 

12.4 

T3 

16.1 

19.6 


TABLE III: Same as Table HI] for combinations of the DO jjdt] and LHCb il8ll data. 


The overall quality of the collider data description achieved in our analysis is sufficiently good, 
cf. Table HVl in Section 11111 However, in cases we observe a certain deterioration of the;^^ values 
as compared to the variants of fit with only one single W-- and Z-boson data set employed in 
combination with the world DIS and the fixed-target DY dimuon data used in the ABM 12 fit. This 
indicates a tension between the data of certain experiments considered. To the most extent, this 
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appears for the LHCb and DO data sets as follows from a comparison of the values obtained 
in variants of our analysis based on including different combinations of the LHCb and DO data of 
Refs. lH, 18, while the other collider data of Table ITVl are omitted, in analogy to the study of 
mutual consistency of the Z-boson LHCb data above (cf. Table HH). The value increases sig - 
nificantly upon combining the DO data on the -boson production [0] with the LHCb ones [fl^ 
(variant #3 of Table Un])^. Additional methodical studies are still necessary in order to resolve the 
tension observed either in favor of the DO or of LHCb experiment. However, in any case, the 
discrepancies are at the level of Icr fluctuations in the value of Therefore the PDF shapes 
preferred by each experiment are basically comparable within the uncertainties. 


III. COMPARISON WITH OTHER DETERMINATIONS 


Constraints on the iso-spin asymmetry in quark distributions coming from the DO Tevatron data 
have also been recently considered in the HERAFitter analysis [Il3n . However, in contrast to our 
case, major parts of the results in Ref. 0 are based on the data for the charge asymmetry of 
VF*-bosons 113 ill , rather than on the lepton asymmetry 113ill . In fact, these two sets stem from the 
same data sample and the VF“-asymmetry is reconstructed from the lepton distributions by solving 
for the kinematics of the VF“-decay with the world-average constraint on the lT“-mass value. This 
solution is not unique and, therefore, the final -distribution is obtained in a statistical way, by 
modeling probabilities of those solutions. The DO modeling is based on the early CTEO PDFs, 
which do not agree with ours in the relevant range of kinematics for the data of Ref. Bill , cf. 
Fig- HI This difference does not allow a consistent use of the DO data on the VF“-asymmetry in our 
analysis. To check this issue in a quantitative way we compare the DO data on the VF“-asymmetry 
with the predictions based on the results of present analysis. Despite a good description of the 
charged lepton data is achieved in our fit, the VF“-asymmetry data systematically undershoot the 
predictions, cf. Fig. HI Therefore the VF“-asymmetry data would introduce a bias to our fit. A 
conceivable impact of those data on the PDFs can roughly be estimated from Fig. It consists 
in a shift of the dju ratio at large x by Icr down and an essential reduction of the error in the 
d/u ratio, since the errors in the data at large pseudo-rapidity are significantly smaller the ones in 
predictions. Such a reduction of errors is also observed in the HERAFitter analysis 0. This 
trend requires further clarification of the DO data statistical significance on the information theory 
framework, since, as we have pointed out above, the same data sample is used to produce both the 
charged-lepton and the VF“-asymmetry distributions. 

Other recent PDF determinations are partially based on the data considered in our analysis, cf. 
Table UVl although a full coverage is not achieved anywhere else ^. Despite this partial coverage 
of data in other PDF fits, the values of;^^ achieved in those fits are generally bigger than ours. In 
particular, this concerns the ATLAS data of Ref. [|^ and the DO electron data of Ref. [01], which 
provide an essential constraints on the quark iso-spin asymmetry in our analysis. The value of;^^ 
obtained in the present analysis for the ATLAS data is also smaller than the one of the ABM 12 
fit d^, evidently due to the relaxed form of the small-x iso-spin asymmetry employed. 

A selected set of the NNLO quark distributions obtained in the present analysis is compared to 
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The problematic data points of LHCb llSil corresponding to the muon pseudo-rapidity tj ^ - 3.375 have been not 
used in the comparisons of Table HIH as we have discussed above. 

The Tevatron data on the charged IT'''-asymmetry are used in the MMHT14 fit ill ill instead. 
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the latest results of the CT NNPDF and MMHT ^ PDF fits in Figs. |9] and [TO] In 
the MMHT 14 analysis the sea iso-spin asymmetry I(x) is parameterized in the Regge-like form 
Eq. ([U), so it vanishes at small x. The CT14 and NN3.0 fits are based on a relaxed form of I{x) at 
small X, and in our analysis. Due to the limited set of the W-- and Z-boson eollider data used in 
those fits, the CT14 and NN3.0 errors in I{x) are significantly bigger than the ones obtained in our 
analysis. However, the CT14 error band for I(x) is comparable with ours, while the one of NN3.0 
goes higher at x ~ 10“^. At x ~ 0.1 all determinations are consistent since this ki nematics range is 
controlled by the fixed-target DY dimuon data, common for all four fits. The shape of I(x) at large 
X in CT14, MMHT14 and our analysis is controlled by the quark sum rule. Therefore it vanishes 
at x = 1. A similar behavior is observed for the NN3.0 result, despite the author’s declaration of 
a model-independent parametrization of PDFs. Evidently, some phase-space constraints are still 
applied in this case. However, a detailed clarification of this issue is impossible in view of the 
fact that the explicit shape of the NN3.0 PDFs is unpublished. The ratio dju at large x obtained 
in our analysis goes lower than other determinations, although it is consistent with them within 
uncertainties, which are much bigger than ours, especially for the MMHT 14 and NN3.0 cases. As 
well as in the case of I{x) this happens since a limited set of the W-- and Z-boson collider data 
is used in those fits. Besides, the cut of > 13 GeV^ on the invariant mass W of the produced 
hadronic system imposed on the DIS data in the CT14, MMHT14, and NN3.0 analyses removes 
the large-x data for the deuteron target, which can provide an additional constraint on the J/m ratio, 
as it happens in the ABM 12 analysis [SI]. Since the d/u ratio at large x is basically uncontrolled in 
the MMHT14 and NN3.0 analyses, its central values differ from ours at x = 0.9 by factor of ~ 30 
and ~ 20, respectively. Such a spread obviously makes these PDFs inapplicable for a precision 
study of physics effects beyond Standard Model at large scales. Meanwhile, the MMHT 14 and 
NN3.0 predictions are in clear disagreement with the forward VF“-boson production data by LHCb 
and DO, cf. Figs. [Hand |3l This means, they can be consolidated with ours once additional W-- 
and Z-boson collider data are included into those analyses, in analogy to the recent CT14 results 
(see the difference between the CTIO and CT14 predictions in Figs. [Hand [3]). 

The techniques used in the other PDF analyses [13-12] to take into account the NNLO correc¬ 
tions to the DY differential distributions are also different. The MMHT 14 analysis is based 

13^ output supplemented by NNLO ^-factors computed with the codes 
33l] . For the most accurate recent data on the lepton asymmetry this ap¬ 
proximation is insufficient since the NNLO corrections demonstrate a dependence on the lepton 
pseudo-rapidity with a spread beyond the data uncertainties, cf. Lig. [^ In the NN3.0 analysis [Il2n 
the NNLO corrections are applied in the form of so-called C-factors Odil including kinematic de¬ 
pendence of the NNLO corrections, however, still computed with the PDLs fixed. Linally, the 
CT14 analysis is based on the ResBos code llisll . which provides the QCD corrections up to 
the NNLL accuracy only. In view of such a variety of tools a consolidation of the PDLs obtained 
in those fits also requires a thorough validation of the accuracy of the theoretical predictions used. 


on the NLO APPLGrid 
FEWZ [El and DYNNLO 


IV. SINGLE-TOP HADRONIC PRODUCTION 


The results for the PDL fit have interesting implications in case of theory predictions for the 
production of single top-quarks at the LHC. The t-channel single (anti)top-quark cross sections cft 


The HERAFitter analysis ill311 is performed to NLO accuracy in QCD only. Therefore it is not considered in this 


comparison. 
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Experiment 

ATLAS 

CMS 

DO 

LHCb 

(TeV) 

7 

7 

1.96 

7 

8 

Final states 

1 T+ ^ tv 

IT“ ^ tv 

Z tt 

1 T+ ^ /i+y 

W- /u-v 

1 T+ 

W- -^/u-v 

1 T+ ^ e+v 

W- e-v 

1 T+ 

W- fi-v 

Z 

Z ^ e+g- 

Reference 


UJ 


L4] 

08] 

[19] 

NDP 

30 

11 

10 

13 

31 

17 


this work 

29.8 

22.5 

16.9 

18.0 

44.1 

18.2 

this work" 

32.3 

19.5(13.5* ) 

13.5 

9.5 

34.7 

19.1 

ABM12 [8] 

34.5 

- 

- 

- 

- 

- 

CT14 [TO] 

_C 

_ d 

- 

34.7 

- 

- 

HERAFitter [13J 

- 

- 

13 

19 

- 

- 

MMHT 14 [llj 

39 

- 

21 

- 

- 

- 

NN3.0 L12J 

35.4 

18.9 

- 

- 

- 

- 


"The variants with all collider DY and VT=‘=-boson data excluded except the one given. 

*The value obtained assuming systematic uncertainties to be uncorrelated. 

"The ATLAS data on and Z-boson production cross sections are used in combination with the data on the lepton 
charge asymmetry. The value of x^INDP - 51/41 = 1.25 is obtained for this sample. 

'^Statistically less significant data with the cut of > 35 GeV are used. 


TABLE IV: The values obtained for the data samples considered in the present analysis, cf. Table U in 
comparison with the ones obtained in other PDF fits. The NNLO results are quoted for all cases, except for 
the HERAFitter results 113]. 


and crj have been measured at the LHC for yfs = 1 TeV by ATLAS iflsll and for yfs = % TeV by 
CMS [[I^ . The inelusive eross seetions at NLO in QCD ean be eonveniently eomputed with the 
Hathor library ISaOTlj, whieh features t-ehannel, 5-ehannel and Wt- single top-quark produetion 
and employs the hard partonie eross seetions at NLO in QCD based on Refs. 113SL13911 for the 
t-ehannel proeess. 

For the PDFs of the present analysis as well as the sets diseussed in the previous Seetion, i.e., 
ABM12, CTIO, CT14, MMHT, NN3.0, we evaluate the eross seetions crt and iXf for the t-ehannel 
produetion of single (anti)top-quarks along wi th th eir ratio Rt = (Xf/crf. The eross seetions are 
eomputed for running top-quark masses in the MS seheme aeeording to the algorithm deseribed 
in Ref. [|4^ . We ehoose mtimt) = 163 GeV in the MS seheme, whieh eorresponds to a pole mass 
^poie _ ]^ 72 . 4 g GeV using the eonversion to the on-shell seheme at three loops, of. Ref. il. 

In Figs, m] and [I2] we display the eross seetions cr^ and cr,- at = 7 TeV and = 8 TeV, 
respeetively, using mt{mt) = 163 GeV at the nominal values for the seales, The 

eentral value for eaeh PDF set is eomplemented by the symmetrized PDF and, if applioable, as un- 
eertainties. It is obvious from Figs. [H] and El that within the eurrent experimental unoertainties all 
predietions agree with data from ATLAS iflSn and CMS Moreover, the results demonstrate 
that the eross seetion predietions are quite stable against higher order perturbative oorreotions. 
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Predictions based on PDF sets taken either at NLO or at NNLO agree with uncertainties. Also, 
the known NNLO QCD corrections [Il4l] to the hard partonic cross section are quite small. Inter¬ 
estingly, in the case with no cut on the transverse momentum of the top quark they are negative, 
with values reported 114] for the ratios -1.6% and -1.3% for pp- 

collisions at = 8 TeV for a pole mass = 173.2 GeV. Note, that the mass dependence of 
the single top-quark cross sections is very mild, cf. Ref. 113711 . so that these conclusions are directly 
applicable to our analysis as well. 

In Fig. [13] we show results for the ratio Rt = crj/cTf of the cross sections. In this quantity, 
theory uncertainties basically cancel, making Rt a very sensitive probe for the dju ratio in PDFs at 
large x. From the ATLAS and CMS data plotted in Fig.[T3]is also evident, that uncertainties due to 
experimental systematics have almost completely canceled, so that current measurements of Rt are 
limited by statistics. All predictions for Rt with the various PDF sets agree with the ATLAS [jlSll 
and CMS ifl^ data within uncertainties, although the results for CTIO, CT14, MMHT and NN3.0 
are shifted systematically to the edge of the Icr interval of the experimental uncertainty band. The 
numbers for ABM 12 and this work do agree very well with the central values for Rt measured by 
ATLAS and CMS. With much higher rates and, as a consequence, improved statistics for top-quark 
production in the high energy run at = 13 TeV, the prospects for single top-quark production 
to serve as a standard candle process are very good. 


V, CONCLUSIONS 

We have used the currently available data for hadro-production of W-- and Z-bosons at the 
LHC and Tevatron from the experiments ATLAS, CMS, DO and LHCb to constrain the light quark 
flavor distributions in the proton. Due to the data’s kinematic range extending to large rapidities 
important new information on the light quark PDFs over a wide range in x has been obtained. 
Given the unprecedented experimental precision with total systematic errors often well below 
1% theory comparisons at NNLO accuracy in QCD are mandatory. The analysis is based on 
fully differential QCD predictions at NNLO for the rapidity distributions of the leptons from the 
W-- and Z-boson decay. The PDF fit requires fast and reliable theory predictions at per mil 
level numerical integration accuracy which poses a challenge for currently available codes. In the 
present analysis the choice has been the publicly available code FEWZ at the expense of typical 
CPU times necessary to accomplish these computations of 0(1 O^h) for each data set. 

The new data considered in the present analysis allow for a refinement of the parameterization 
for the u- and J-quark distributions. The PDF fit has been performed as a variant of the ABM12 
analysis, which is based on the world DIS data, on DY dimuon data from fixed targets and on early 
LHC data for hadro-production of W-- and Z-bosons. Accommodation of the new data in the fit 
can be achieved by exploiting the parametric freedom either in J-quark distribution at large x or in 
the J-quark distribution at small x. The agreement with data for forward -boson production at 
large rapidities rj ~ 2.5 3 for DO and in the range 77 ~ 3 4.5 for LHCb has suggested a modifi¬ 
cation of the J-quark distribution at small x. By relaxing the commonly chosen parametric ansatz 
for the iso-spin asymmetric sea x{d - u) at small x, both DO and LHCb data can be consistently 
accommodated in the fit. This leads to a non-vanishing x(d-u) asymmetry at values x ~ 10“"^ and 
a delayed onset of the generally anticipated Regge asymptotics at small x, a feature which is not 
entirely without precedence in perturbative QCD. 

In the analysis great care has been taken to consider only data sets which are mutually consistent 
and we have demonstrated that there are issues with the LHCb data. In checking the consistency 
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of LHCb data for Z-boson production, we have noted tensions among the various data sets, whieh 
has lead to the exelusion of the Z ^ e^e~ LHCb data taken at = 7 TeV from the fit. We 
have also observed that different data for the deeay of forward -bosons published by LHCb 
are not eonsistent with eaeh other regarding the shape and the treatment of final state radiation 
to aeeount for QED eorreetions. This has lead us to reeonsider the quoted uneertainties for this 
souree. Finally, the LHCb data for VL“-boson produetion display some tension with the respeetive 
DO data sets, whieh is notieeable from a deterioration of the values in the eombined fit. In 
that ease, however, eurrent information is insuffieient to resolve the observed diserepaneies whieh 
show up as fluetuations at the level of Icr in favor either of the DO or the LHCb experiment. 

The results of the present analysis and the data have been eompared to predietions obtained with 
other global PDF sets. The theory predietions at NNFO in QCD for the eharged lepton asymmetry 
from VT“-deoays based on CTIO, MMHT or the so-ealled unbiased fit NN3.0 undershoot the data 
at large rapidities dramatieally, leading to poor eompatibility of those predietions with DO and 
FHCb data in the forward region. This indieates a signifieant eonstraining power of the data for 
those global fits, as has been pointed out in a reeent HERAFitter analysis as well. 

As an applieation to eollider phenomenology, the fit results have been used to prediet the eross 
seetions for the produetion of single (anti)top-quarks at the FHC, cr,, aj and their ratio Rt = CtlcTh 
demonstrating very good agreement with existing data and underpinning the potential of this pro- 
eess to serve as a standard eandle. 

In summary, the analysis shows the gain in knowledge on the light flavor PDFs from eurrent 
FHC and Tevatron data. We expeet further improvements from data on hadro-produetion of W-- 
and Z-bosons as well as from high statisties measurements of single (anti)top-quark produetion 
in run II of the FHC at = 13 TeV. Additional advanees in the small-x regime ean be expeeted 
from the proposed Farge Hadron Fleetron Collider (FHeC) [Id^l . whieh is to eollide an eleetron 
beam up to possibly 140 GeV with the intense hadron beams of the FHC, ineluding deuterons, and 
whieh extends the eurrently eovered kinematie range in Bjorken-x signifieantly. 
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DO (1.96 TeV) 



FIG. 1: The pulls of the new DO data on the charged lepton asymmetry for pp —> W- + X —> for elec¬ 
trons 10] (circles) and muons |0] (squares) at the Tevatron with yfs = 1.96 TeV as a function of the pseudo¬ 
lepton rapidity rji normalized to the predictions with FEWZ (version 3.1) at NNLO in QCD computed in 
the present analysis. For comparison, the results of the ABM12 fit with the corresponding Icr band of the 
combined PDF-t-a^ uncertainty and the central fits of CTIO CT14, MMHT, and NN3.0 are displayed. 
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LHCb (7 TeV, 1 fb 



FIG. 2: Same as Fig. [T] for the LHCb data 1181] on the production of (left) and lF“-bosons (right) in 
the pp collision at ^fs = 1 TeV (circles) in comparison with those data before a trigger efficiency has been 
fixed 10] (squares). 
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FIG. 3: Same as Fig.[T]for the new CMS 111] and LHCb data 0] on the muon charge asymmetry in inclusive 
pp —> W- + X ^ p-y production at the LHC with sfs = 1 TeV. 
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LHCb 



FIG. 4: Same as Fig. [T] for the LHCb data on the Z-boson production il8Lll9L 12311 at the c.m.s. energy of 


7 TeV (left) and 8 TeV (right) in the muon channel (circles) and the electron channel (squares) in comparison 
to the ABM12 predictions. 


LHCb 



FIG. 5: Left panel: Correction for the final-state-radiation (FSR) effects applied in the LHCb analysis lll8h 
to the data on the pseudo-rapidity distribution of muons produced in forward VF*-decays (dashes: , dots: 
jT, solid: difference befween fhese fwo). Righf panel: The same for fhe earlier LHCb analysis of Ref. 12ih . 
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D0(1.96 TeV, 9.7 fb 
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FIG. 6: The experimental uncertainties in the DO data on the electron asymmetry |0] (points) compared 
with the numerical integration accuracy achieved with FEWZ (version 3.1) 11241. 12511 in the generation of the 
NNLO grids employed in present analysis (shaded area). The impact of the NNLO corrections on the hard 
scattering coefficient functions is indicated by the dashed curve. 


N^4, |x=3 GeV 



X X 



FIG. 7: Left: The iso-spin asymmetry of the sea I{x) in the 4-flavor scheme at the factorization scale 
/I = 3 GeV obtained in the present variant of the ABM12 analysis (solid lines: central value, dots: Icr error 
band) in comparison to the Icr band of l{x) obtained in the ABM 12 fit fsl] (shaded area). Right: Same for 
the ratio dju at large x. 
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FIG. 8: Same as Fig. [T]for the DO data on the charged IF*-asymmetry Bill extracted from the DO data on 
the electron charge asymmetry 0] as a function of the IF*-rapidity rjyv. The shaded area displays the PDF 
uncertainties in the predictions. 
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N^4, ^=2 GeV 



FIG. 9: The Icr band for the iso-spin asymmetry of the sea I(x) for the 4-flavor scheme at the scale of 
ju = 2 GeV as a function of the Bjorken v obtained in the present fit (gray shaded area)in comparison 
with the correspon ding ones obtained in the CT14 lE^ (red right-tilted hatch), MMHT14 IJT] (blue dashed 
lines), and NN3.0 il2l] (green left-tilted hatch) analyses. 
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FIG. 10: Same as in Fig.|9]for the ratio dju. 
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FIG. 11: Cross sections together with their Icr PDF uncertainties for the t-channel production of sin¬ 
gle (anti)top-quarks in pp collision at = 7 TeV in comparison to ATLAS data ilSll for a MS mass 
nttimt) = 163 GeV at the scale pR=pf = mt{mt) with PDF sets are taken at NNLO. The inner (yellow) band 
denotes the statistical uncertainty and the outer (green) band the combined uncertainty due to statistics and 
systematics. 
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FIG. 12: Same as Fig.[TT]for pp collision at = 8 TeV in comparison to CMS data I16ll . 



FIG. 13: Same as in Fig.[IT]for the ratio of cross sections Rt = (Ttlcrj in comparison to ATLAS data ilSll at 
yfs -1 TeV (left) and to CMS data |16] at ^[s = % TeV (right). 
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